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1. Introduction 

Since the pioneer works of Jost and Pais in 1951 pQ and Bargmann in 1952 
the determination of upper and lower limits on the number of bound states of a 
given potential, having spherical symmetry V(r), in the framework of non relativistic 
quantum mechanics is still of interest. A fairly large number of results of this 
kind can be found in the literature for the Schrodinger equation (see for example 
[31 El El E3 El El El EH [HI H21 CHI HH HSl EH C3 US! and for results applicable to 
one and two dimension spaces see for example 1201 EH 1221 EH] ) 

An important theorem for classifying these results was found by Chadan |S| and 
gives the asymptotic behavior of the number of £-wave bound states as the strength, g, 
of the central potential V(r) = gv{r) goes to infinity: 

„l/2 „oo 

NtFH- — / drv~(r) 1/2 as g -> oo, (1) 

7T JO 

where the symbol « means asymptotic equality, V~(r) = gv~(r) and v~{r) = 
max(0, —v(r)) (see also Ref. [21] for a generalization of relation (JTJ). This result implies 
that any upper and lower limit on N( which could yield cogent results should behave 
asymptotically as g 1 ^ 2 . More importantly, the relation (JTJ gives the functional of the 
potential, that is to say, the coefficient in front of g 1//2 that appears in the asymptotic 
behavior. Upper and lower limits on the number of £-wave bound states featuring the 
correct g 1 ^ 2 dependency was first obtained in Ref. [Hj- Upper and lower limits on 
N( featuring the correct asymptotic behavior (JTJ) was first derived in Refs. [T71 HHj . 
In practice, the asymptotic regime is reached very quickly when the strength of the 
potential is large enough to bind two or three bound states. 

The situation is completely different when one consider the transition between zero 
and one bound state. In contrast to one and two dimension cases where any attractive 
potential, satisfying adequate integrability conditions, has at least one bound state, in 
the three dimensional case, the potential acquires a bound state only if it is attractive 
(negative) enough. Thus, for central potential for example, there exists a "critical" 
value, gg\ of the coupling constant (strength), g, of the potential, V(r) = gv(r), for 
which a first £-wave bound state appears. The determination of this critical value 
requires to solve the zero energy Schrodinger equation [31 1213 EE]- To circumvent the 
exact calculation of the Jost function at zero energy, upper and lower bounds are very 
helpful. From now on, use will be made of the standard quantum-mechanical units 
h — 2/i = 1 where /x is the reduced mass of the particles. 

In 1976, Glaser et al. have obtained a strong necessary condition for the existence 
of a £-wave bound state in an arbitrary central potential in three dimensions |9 j : 

This relation yields a lower limit on the critical value gff), by making a minimization 
over p > 1, which was shown to be very accurate (see for example llB*! I28j). 
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Other necessary conditions for the existence of bound states can be found in the 
literature (see for example [7HJTJI2E] and for reviews see [JUJ HZ1 HE] ) , but in general, the 
relation (J2J) yields the strongest restriction on gf® (in some cases, the relations obtained 
in Ref. [2Hj can however be better). 

Sufficient conditions for existence of bound states, yielding upper limits on the 
critical value of the strength of the potential, are scarcer. Two sufficient conditions for 
the existence of at least one bound state with angular momentum i have been found by 
Calogero in 1965 jSHEj 

pR poo 

(VR>0) / drr\V(r)\(r/R) 2e+1 + / dr r \V(r)\ (r I Ry {2l+1) > 2£ + 1, (3) 

JO JR 

and 

poo r -, _1 

(Vi?>0) R dr\V(r)\\(r/R) 2t + (r/R)- 2e R 2 \V(r)\\ > 1. (4) 
Jo 1 J 

These two conditions apply provided the potential is nowhere positive, V(r) = — \V(r)\. 

The most stringent conditions obtain by minimizing the left-hand sides of (jHJ) and 

over all positive values of R. Some other sufficient conditions for the existence of bound 

states can been found in the literature (see for example [TBI EH] and for reviews see 

^OJEIE])- A sufficient condition which does not require the spherical symmetry for 

the potential V was proposed in 1980 by Chadan [SH] (see also jSHl)- When the potential 

is central and purely attractive, the inequality: 

TtK^>TtK^\ (5) 

where 

K {1 \r,r') = inf (r,r') m sup(r,r')"* V(r') 

POO 

K (n \r,r')= ds K {1 \r,s) {s,r'), (6) 

J 

implies the existence of at least a bound state for the potential V(r). In the case where 
V(r) has some changes of sign, the condition (J3J, is replaced by 

TrK (4) >M (2) , (7) 

which implies that one of the potentials ±V(r) has at least one bound state. 

Recently, an upper limit on the critical strength has been found, originating from 
a variational technique |31j : 



(8) 



poo r poo py 

gW < A / dxF(2p-l;x) / dyF(p;y)y- x dzF(p;z)z 
Jo Uo Jo 

with F(q;x) = x q v(xp q+1 '/ 2 , v(x) > 0, A = i + 1/2 and q > which was found to 
be very accurate. Clearly more accurate upper limits could be obtained but depending 
strongly on the choice of the trial wave function. 

The present article follows a different scheme. To circumvent the difficulty to 
guess a trial (strongly potential dependent) wave function for variational methods 
we propose upper and lower limits originating from iterative procedures, all designed 
to converge towards the exact result. The methods proposed by Chadan enter this 
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category of iterative convergent procedures. Also, Lassaut and Lombard [16^ worked 
in this sense some years ago, but the procedure was constrained by the condition 
gf> J °° drr\v(r)\ < 2(2£ + 1), which is not verified when the convexity of the potential 
is too high or when the angular momentum I increases. We thus propose in this article 
sequences of upper and lower limits of the critical value, gf\ which converge towards 
the critical coupling constant without any restriction on the possible values of gf>. The 
advantage of our procedure, based upon the Riesz theorem [321 IHH] m what concerns the 
upper limits, is that there is no need to start from a conveniently chosen wave function. 
Indeed, the improvements are simple for these sequences: one just need to calculate the 
next order. Note that the basic idea of sequences of lower limits for have already 
been explored in Ref. [2%] . 

The paper is organized as follows. In section |2] we derive the upper limits on the 
critical value gjfK In section 0] the algorithm for generating both lower and upper limits 
is discussed. In section |U our proposal of upper and lower bounds are tested against the 
exact values for common potentials. Our conclusions are presented in section El 



2. Upper limits on the critical strength 

From now on, we assume that V(r) is locally integrable and such that 

POO 

/ drr \V(r) I < oo, (9) 
Jo 

remembering that we consider purely attractive potentials namely satisfying V(r) < 0. 
Following Birman and Schwinger 123 the critical values of the strength of the 
potential correspond to the occurrence of an eigenstate with a vanishing energy. In this 
paper we consider the zero energy Schrodinger equation that we write into the form of 
an integral equation incorporating the boundary conditions 

POO 

ue(r) = - dr' ge(r,r')V(r')ut(r'), (10) 



o 



where gi(r, r') is the Green function of the kinetic energy operator and is explicitly given 
by 

9t(ry) = ^r{^ r~\ (11) 

where r< = inf (r, r') and r> = sup (r, r'). An important technical difficulty appears if 
the potential possesses some changes of sign (see relation ()12j) below). This is overcome 
for the derivation of necessary conditions, or of upper bounds on the number of bound 
states, by replacing the potential by its negative part V(r) — > V~(r) = max(0, — V(r)). 
Indeed, the potential V~(r) is more attractive than V(r) and thus a necessary condition 
for existence of bound states in V~(r) is certainly a valid necessary condition for V(r). 
This procedure can no longer be used to obtain sufficient conditions. For this reason we 
consider potentials that are nowhere positive, V(r) = —gv(r), with v(r) > 0. 

To obtain a symmetrical kernel we now introduce the function ipi(r) as follow 

Mr) = ^fiirjuiir). (12) 
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Equation (JTUJ) becomes 

roo 

ipe(r) = g dr' K e (r,r')ipe(r'), (13) 
Jo 

where the symmetric kernel Ki(r,r') is given by 

K t {r, r') = ^fi{r)gt{r, r') y^V). (14) 

The relation (fT3*|) is thus an eigenvalue problem with a symmetric kernel and, for each 
value of £, the smallest characteristic numbers are just the critical values g^>. The 
higher characteristic numbers correspond to the critical values of the strength for which 
a second, a third, £-wave bound state appears. The kernel (fl~4"|) acting on the Hilbert 
space L 2 (M) is an Hilbert-Schmidt kernel [31] for the class of potentials defined by (JHJ), 
i.e. satisfies the inequality 

roo roo 

/ / dxdy Kt(x,y)KAx,y) < oo. (15) 
Jo Jo 

Consequently the eigenvalue problem (fTB^l always possesses at least one characteristic 
number [35] (in general, this problem has an infinity of characteristic numbers). 
We propose now to solve the eigenvalue problem (|T3|) using iterative methods. 



2.1. Iterative power method 

Let us write, for simplicity, the relation (fT3*|) under the form 

ip e = gJCeipe, (16) 

where /Q denotes the symmetric linear operator, operating on the Hilbert space L 2 (M), 
which is in this paper the integral operator generated by the so-called Birman-Schwinger 
HUES] kernel K t , Eq. (JUJ). 

Since the kernel Kg(r, r') is Hilbert-Schmidt, /Q is a compact operator |27j . As /Q is 
symmetric the Riesz theorem applies |32] 133] . For each value of the angular momentum 
£, the set of eigenvalues l/g p , 1 < p, (which in the present case are all positive) can be 
ordered according to a sequence tending to zero, 1/gi > l/g 2 > ... > l/g p > ... > 0. 
There exists an orthonormal basis in L 2 (IR), labeled v 9 p( r )> V ^ 1> each <p p (r) being 
associated to l/g p , and for each function 0^(r) e L 2 (M) 

oo oo oo ^ 

p=i p=i p=i 9p 

where the symbol (f\g) denotes the scalar product J °° dr f(r) g(r). For the sake of 
simplicity we have dropped the indices (£) which should appear on g p and on ip p (r). 
The positivity of the eigenvalues originates [33J from the fact that /Q is positive i.e. 

(V0,GL 2 (M)) (&|£/&>> 0. (18) 

In the present case, where the potential has the spherical symmetry, there is no 
degeneracy and we have strict inequalities for the eigenvalues 
111 

— > — >...> — > . . . > 0. (19) 

9i 92 g P 
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This is due to the fact that the eigenstates are solutions of a linear second order equation 
with constraints at the origin and infinity. 

Now, we introduce the iterated kernel Kg(s,t) of Kf(s,t) 

Ki n \s,t) = r ' duKAs,u)K [ I l ~ l \u,t) n>2, (20) 
J o 

with 

K^(s,t) = K £ (s,t). (21) 
We can then compute the scalar product between <p£ and <fie, and find 

OO -1 

(&|XtV<> = <4^I<M = £ - (0,|^) 2 . (22) 
Therefore, we obtain the following convexity-type relation 

2 



p< g =l \9p9q) \9p 9q, 



X 



> 0. (23) 
Since the left hand side of this latter equality is positive we conclude that the sequence 

l *t +1) ^) 

i*tv<> 



with 



6o = (25) 



is always increasing. 

On the other hand, taking into account (}2*2*|) and (J23jl we have 



9i E P °°=i(S) n ^l^) 



Due to the strict inequalities ()19|) and the Parseval's formula 

OO 

11^112 = Etota) 2 . (27) 

p=l 

the sequence <5 n converges to l/<7i, where g\ = g^ is the lowest critical value of the 
coupling constant of the potential, except when {4>i\ipi) is zero. This procedure, known 
in the literature as the iterated power method, yields the maximal eigenvalue of the 
problem considered, i.e. in the present case to the lowest critical value gjf) = g\. Any 
starting positive (non zero) squared integrable function (f>i(r) is appropriated. Indeed 
the function </?i(r) has no node which insures that the scalar product (0^|<^i) is not 
zero. (In our numerical studies, presented in section 0] use is made of the choice 
4>e(r) = r l+1 ^Jv{r)). Moreover since n t— > 5 n is increasing we get upper bounds for 
gf> namely gf> < ... < 1/8 P < ... < l/5 2 < l/8\. 

There exists other iterative methods in the literature and we discuss briefly two 
variants of them in the next two sections. 
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2.2. Kellogg 's method 

In this section, we consider the method proposed by Kellogg for the compact operator 
K.£ jHEl E3] • We construct the following sequence of functions 

4> < f +1 \r)= I" dr , K t (r,r r )^ ) (r / ), (28) 
Jo 

where K^(r,r') is given by ()14j) . This latter relation is schematically written as: 

0f +1 ) = /Q0f. (29a) 
For non zero (f>i(r), Kellogg considers the following sequence of numbers 

110. 



ln + l= n >!> ( 29& ) 



where ||<^ n ' ) ||2 is the L 2 norm of the function tp^ir). 

We keep our conventions namely y?i(r), f2{f)i still denote the eigenfunctions of 
the problem (jl6j) and g\ < g% < . . . are the corresponding characteristic numbers. 

Suppose that 4>f\r) is orthogonal to the functions v?i( r ); V 9 2( r ), V 9 fc-i( r ) but not 
orthogonal to the function (pk( r )- Then the sequence 7„ converges toward g% with the 
property that g% < 7„. Moreover the sequence of functions (j)^^) / ll^^ \h converges 
to v 3 fc( r ) m £ 2 (R) [HS]. Consequently there exists a subsequence of n 1— > ^ (r)/||<^ || 2 
which converges almost everywhere to (pk(r). 

The convergence of 7„ is illustrated simply here, where we are interested by the 
smallest characteristic number g^K Still we choose a positive (non zero) squared 
integrable function <pf\r), which, according to our previous discussion, is not orthogonal 
to fi{r). The sequence of numbers 7„ provides us upper bounds for because this 
sequence is decreasing and converges to gj®. The monotony of the sequence n 1— > 7„ is 
related to the following equation, derived from (|22|). 

ll^ll 2 . = (^Uf\4 n ¥" ] ) = E 4r <*f I**) 2 . (30) 

p=i yp 

which leads to the inequality 

.,2 _„,2 MU(«+2)||2 iu(n+l)||a_ _ 1 (jL_M /^(0)|,„ \2/^(0)i 



(7^2-7^)11^11^ imt'iih -^(^^-^J <CKm> 9 > 2 , 

< . (31) 

The positivity of 7 n asserts 7 n+2 < 7n+i- 

On the other hand we can check easily that 7„ converges to gjp for n going to 
infinity. Indeed, taking into account definitions (j29fc|) and f|3(J|) we have 



sr=i(g) 2n (0fl^) 2 



2 _ 2 ^P=^g„' ^ /OON 



if 

which, using (|T9*|) and (|2*7|). shows clearly the convergence. 
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In the section^ use will be made of the iterative procedure (J28|) but for the functions 
u^\r), defined by u^\r) = <j>^\r)/ 'Jv(r), namely 

u*P +l \r) = / dr' v(r')gi(r, r')«| (r'), 



dr' r' t+1 v{r')uf\r') + H dr' r'- l v(r')uf\r'), (33) 



with (r) = r e+1 . The sequence 7 n is given by 



7n+l 



/ °°rfrt»(r) [4 n) (r)] 2 ^ 



\ f^drv(r) [uf +1) (r)f 



2.3. Kolomy's method 

The second variant of iterative methods was proposed by Kolomy jHZj who constructed 
the sequence of functions 

,(n+i) _ {K>i4>t \4>t ) r An) 



1 Mr ^a(^) i 1 2 

W^lVl \\2 



■/Q0f j , (35a) 



and consider the sequence of numbers 

The factor ((/Q</>^|</>^))/(||/Q0^|||) in (j35aj) is useless as far as the series of number 
(1356)1 is concerned. In this paper, where our interest is focused to the lowest eigenvalue 
of the problem, we simply drop the normalization factor in (J35a)l , namely we consider 

# +1) = ^# } , (36) 

coupled with ()356j) . 

When the starting trial function (j)f\r) is a positive (non zero) squared integrable 
function, it is not orthogonal to (pi(r). Then the sequence f3 n converge to g\ with 
9i < Pn- We still obtain a sequence of upper bounds for gf\ Indeed the sequence 
n i— > /? n is decreasing since 

(^-/WOII^^II! = - ^(^r t-^) 2 + >< 

< 0. (37) 

On the other hand we can check easily that (3 n converges to g\ = g$p for n going 
to infinity. Indeed, taking into account definition (j356j) and (j!7j). we have 

which, using (fE?)) and ()27)1. shows clearly the convergence. 
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For the sake of convenience, we rewrite the iterative procedure (|36|) in terms of the 
function (r), defined by u^\r) = (fr t n ' (r)/Jv(r). Using (|356jl . the sequence (3 n reads 



Io° drv{r)uf\r)uf +1 \r) 

f£°drv(r) ['U^™ +1 ^(r)] 2 

As we show in section IH both iterative procedures depicted in 12.21 and 12.31 converge 
very rapidly, when use is made of the initial function (pf\r) = r e+1 ^Jv(r). However to 

decrease the number of iterations, more flexible function (ffi (r) , obeying to the desired 
conditions, positivity and squared integrability, could be very helpful. This is explored 
in the next section. 

2.4- Combination of iterative and variational methods 

The variational method we propose is based upon the theorem EH] which states that 
for a symmetric compact operator, 

m V [(JC t W)] = -± 1 , (40) 
i> 1 01 1 

under the constraint that ip is squared integrable and normalized to unity (HV'lh = !)■ 
Since in this article the kernel we consider is positive we have \gi\ = g\. The maximal 
value of (jlOJ) is reached for ip{r) = <fi(r), the eigenfunction associated to the smallest 
eigenvalue g 1 = g^'. Consequently, for any function ip( r ), normalized to unity, we obtain 
the following upper limit 

gPzVbm- 1 - (41) 

This method was used in Ref. [3T] to obtain the upper limit (jHJ) using the trial 
function 

i)(r) = A [r 2p -yr) p ] 1/2 , p > 0. (42) 

The degree of flexibility increases with 

ij){r) = A[r p v(r) q } 1/2 , (43) 

where still (p,q) are varied in the domain ensuring that ip(r) e L 2 (R) and A is a 
normalization factor designed to have \\1pW2 = 1- When we combine the iterative 
methods together with the variational one, we construct a sequence of functions which 
automatically converges in L 2 (M) to the exact zero energy wave function. This means 
that the variational method PD*|) can be used with all iterated ip^ = )C^ip, i = 1,2,..., 
of V(r), Eq. ipjl. 

These possibilities are tested in section 0] and shown to greatly improve the 
convergence towards the exact result. 
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Now, let us introduce another method which provide both lower and upper limits of the 
critical value required. Note that this method has some link with the iterative power 
method discussed above. 

In this section we restrict to the S-wave case since the determination of the critical 
value for the potential V(r) in the £-wave is equivalent to the determination of the 
critical value for the potential 

1 V ( r i/(2e+i)\ 

WAr) = — - (44) 

vy ^'> (21+ 1)2 rttW+i) y ' 

in the S-wave [To] . 

Theorem Let V(r) be a central potential with V(r) = —gv(r),v(r) > 0. Let the 
functions ip n (r) defined by the following recurrence relation 

POO 

"0n( r ) = / dr' g(r,r')v(r') 4> n -i( r ') n>\ (45) 



with ipo(r) = r and g(r,r') = inf(r, r'). Let a n and uo n be two sequences defined as 
follows, 

a n = lim ^ w+l(r) = fn+M = ^ drv ^^(r) n 
r^o il) n (r) i// n (0) f™drv(r) ip n -i(r) 

where the prime denotes the derivative with respect to r and 

Un = nm ^ + i{r) = JS°drrv(r)Mr) n > 1 (4?) 
r^oo ^ n ( r ) / oc drr«(r)V'n-i(»") 

Then, the two sequences a n and cu n converges to l/g(°\ g^> being the critical value of 
the strength of the potential V(r) for which a first S-wave bound state appears. The a n 
sequence is decreasing while the u n sequence is increasing. 

The proof that u n is increasing and converges towards the desired value is rather 
simple. The change of function 0o^( r ) = \Jv{r) ip n (r) in the iterative process (|4*H|) leads 
to 

0<j n+1) (r) = r dr'K (r, r') <^ n) (r') n > 1, (48) 
Jo 

with <f) °\r) = ry / t>(r) and Ko(r,r') = yjv(r) inf(r, r')Jv(r'). Clearly we recover 

# +1) =^o^ n) = 4 n+1) 4 0) n>0, (49) 
in terms of the iterated kernel K.^ (see section E3} - As 

roo 

lim ip n {r) = / drr v(r) ip n -i(r) n > 1, (50) 
is in fact the scalar product 

hm^(r) = (#|^- 1) ), (51) 



we deduce that tu n = 5 n _i when the starting function in (j24j) is 0o( r ) = r^Jv(r). This 
shows that tu n is an increasing sequence converging to l/g^- 
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Now, let us consider the sequence a n . We show in the Appendix A that 



n i — ► ot r 



is a decreasing sequence which converges to l/g^°\ In the |Appcndix A we relate all the 



entering the definition of a n to the Jost function for V(r) at zero energy. This is 
of some interest in the measure where we obtain lower limits for g^ circumventing the 
constraint / °° dr r\v(r)\ < 2 of [Tfij . 

The sequence a n and u) n are simple enough to write the first members explicitly. 
First of all we remark that 

POO 

cto = dxxv(x), (52) 
Jo 

which is just the Bargmann-Schwinger necessary condition for the existence of bound 
states. We also have u> = 0. The next order yields the following relation 

J^dxxv(x) 1 (o)l J °°dxx 2 v(x) 

— 9c — — o poo i / x rx i 5 TT' \ 06 ) 



f™dxv(x)f™dy M(x,y)yv(y) a x c u x 2 / °° dx xv(x) Jq dyy 2 v(y)' 
4. Tests of the bounds on some common potentials 

In this section, we propose to test the accuracy of the various upper and lower limits 
obtained in the sections 121 and El with three potentials: a square well potential (denotes 
hereafter SW potential), 

V(r) = -gR- 2 6(l -r/R); (54) 
an exponential potential (denotes hereafter E potential) 

V(r) = -gR~ 2 exp(-r/ J R); (55) 
a non monotonic potential (denotes hereafter PE potential) 

V(r) = -^ 3 rexp(-r/ J R); (56) 

In these potentials, the radius R is arbitrary (but positive). Due to the scaling 
property, which does not affect the critical value gg\ the radius R appearing in the 
potential can be set to unity. 

In order to test the reliability of the different methods when the convexity of the 
potential increases, we perform simple analytical calculations (that we do not report 
here) involving the SW potential (J54j) . when the angular momentum i increases. We 
also study the E and the PE potentials only for the S-wave. 

In order to examine in what extent our recursive procedure lead rapidly to the exact 
result, the exact value of the critical coupling constants of the potentials, gf \ and its 
approximated upper and lower limits investigated in the sections El and 01 are depicted in 
the Tables IaTI to IA4I More precisely, the Kellogg's coefficients 7 n are shown in the Table 
IA1I the Kolomy's coefficient are reported in the Table IA2I The results obtained with 
the variational method and its iterated are given in Table IA3I In Table IA4| we report 
the results obtained with lower and upper limits discussed in section El As a further 
information, we give in Table lA"5l the results for and uj^ which are compared to those 
obtained with the formulas (J2J), © and (jHJ) for the three potentials (|H4*j) - (|5l)j) . 

In all cases, few iterations are enough to obtain strong restrictions on the possible 
values of g^\ especially for low value of the angular momentum. 
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5. Conclusions 



In this article we have presented in section El several iterative procedures yielding 
sequences of upper limits on the critical value, gf\ of the coupling constant (strength), 
g, of the potential, V(r) = —gv(r), for which a first £-wave bound state appears. In 
section 01 we have obtained a method yielding both upper and lower limits on g^> . All 
these sequences converge rather rapidly to the exact critical value as shown by the tests 
presented in sectionH] Due to the construction of the sequences, for example (3 n = l/^n, 
the convergence of the sequences obtained in section |21 is faster than the ones of section 
01 However results which are the closest to the exact values are obtained, with minimal 
numerical efforts, from the combination of variational and iterative methods. 

Note that the results presented in this article, of similar accuracy than the results 
given by other methods, can always be improved when use is made of a supplementary 
iteration. The accuracy of the results that we obtain can be measured by the difference 
l/u n - l/a n . 



Appendix A. Monotony of the sequence a n 

In this appendix we study the sequence a n defined by pfijl. As made for uj n , we express 
a n in terms of a scalar product. We know that ip n+ i{r) ~ r / °° drv(r) ip n (r),n > 1 at 
the vicinity of r = 0, so that its first derivative is the scalar product 

lCi(O) = (Qo\^) (A.l) 



with qo(r) = \Jv(r) and 0o^(r) = ^Jv(r)ijj n (r). 

The proof of the monotony is made in two steps. First we show that the inverse of 
the Jost function at zero energy, given for example in Eqs. (13) and (14) of Ref. |lfij . 
in terms of the coupling constant g is simply the series 

1 oo 

tsmtS*^ (a - 2) 

which converges for g < g(°>. It can be easily verified by considering ip n (r), defined 
by the recurrence relation (J45|) but where the first term ipo(r) is equal to unity. The 
equation becomes 



tp n (r) = M n — J dr' (r' — r) v(r) ^ n _x(r') n > 1 

M n = / drrv(r)ip n -.i(r) n > 1. (A. 3) 



Introducing 0o^( r ) = $n(r) \] v ( r )i ( t ) o l \ r ) = V , n( r ) \j v {. r ) an d = Jv{r) we 

have, due to the symmetry of the scalar product, 

= (qo\4 n ~ 1] ) = rdrv{r)Mi{r) = <(0). (A.4) 



This allows to define M„ for n = and we have Mn = 1. 
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On the other hand, it has been shown that the zero energy Jost function could be 
written as [TH] 

oo 

/oMHERX^"' (A.5) 

n=0 

with 

/■OO /"OO /"OO 

a n = dr x 7*1 w(ri) / dr 2 (r 2 - ri)v(r 2 ) . . . / dr n (r n - r n ^)v(r n ) n>2 (A. 6) 

JO Jr\ Jr n —i 

and CLq — 1, OL\ = J °° cfr ru(r). We have the following relation between M n and a n 

n 

J]M n _ p a p (-) p = n>l, (A.7) 

whereas for n = 0, ciqMo = 1- This shows that the series J2%Lo M n g n is equal to l//(0, g>) 
for g < Indeed, the relation (jA.7|) is just the relation that exists between Taylor's 
coefficients of the functions f(x) and 1/ f(x). 

Secondly we prove the following lemma 

lemma Let be the series 

oo 

F(g) = Y,Png n (A.8) 

n=0 

convergent for g < R, (R > 0), such that (3q = 1 and 
(Vn > 0) (3 n > 

(Vn>0) P n+2 Pn<P 2 n+1 . (A.9) 
Then, for every A > ; the series 

oo oo oo 

G(<7) = Y,Pp9 p x^AV = E7n^ n (A. 10) 

p=0 g=0 n=0 

convergent for g < R = inf(i?, 1/A) is such that 7o = 1 and has the property 
(Vn > 0) 7n > 

(V72>0) 7n+2 7n<7n+l- (A.H) 

The proof of the lemma originates from the definition of 7„ 

n 

7n = £A^ n " P (A- 12) 

p=0 

which asserts that 70 = 1 and 7„ > since A and /?„ are positive (see (|A.9|l ). The second 
property in (jA.llJ) is satisfied if and only if the following relation is satisfied for n > 

n+2 n /n+1 \ 2 

£ & A n+2 ~ p £ 0, A™ -9 - [J2 P P A n+1 - p < 0. (A. 13) 

p=0 g=0 \P=0 / 

The inequality ()A.13|) is equivalent to the requirement that 

n 

(P n+2 - (3 n+1 \) P P A n " p - Pl+i < 0, (A.14) 

p=0 
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which is manifestly satisfied when 

(Vn, 1 < p < n) (3 n+2 (3 p - (3 n+1 (3 p+l < 0. (A. 15) 

The relation (|A.9J) and the positivity of (3 n imply that this latter inequality, which 
originates from the iterated of (|A.9|) when n is lowered up to p, is verified. 

Now we use the fact that, for any g < g x = gjP> the Jost function can be written as 

m 

fo(9, o) = ( i - f) (A-16) 

n=l \ 9nJ 

where the g n y s still denote the characteristic numbers of the eigenvalue problem 
considered. The product (jA.16|) exists for g < g\ when each series Y^Lo 9n P converges 
for every integer p > 1, which is true when the (positive) trace of the iterated kernel 
Kq (r, r') is finite. The equation (jA.16|) implies that for any g < gi we have 



lim S 



N 



fo(9,Q) n=lp^0\9nj 

Sn =nzf ■ ( A - 17 ) 

n=lp=0 \yn/ 

Since all the quantities in (jA.17|) are positive we can write, for g < g^. 

oo 

Sn = T,4 N) 9 p (A. 18) 

p=0 

where Sq = 1. Now assuming that for some N, the following property holds: 

(Vp>0) 4^ - (4+i) 2 < (A-19) 
according to the lemma, the property still holds for iV + 1 as well. Note that the radius 
of convergence R of the lemma is always minorated by g\ > 0. For N — 1, is simply 
l/g p and the property (jA.19|) is valid. For N = 2, = Y^k=o 9\ k 92~ P anc i ^ can be 
verified that the property (|A.19|) is again valid. Therefore for every N > 2 the property 
(IA.19|) is also valid and in particular for iV going to infinity. The comparison between 
(TOJ) and (IA~T7|) shows that ^(0) = • From the relation (IA~T9|) we obtain 

(Vn > 0) < %#> (A.20) 

which, with the definition (|4T)j) of a n , prove that the sequence a n is decreasing. 

Since a n is decreasing and positive it converges towards some a > 0. On 
the other hand we know that the radius of convergence of the series l/f(g,0) is 
g\. Using the d'Alembert rule for the series of positive numbers (^(0))„>o we have 
1/gi = lim^oo ^ +1 (0)/^(0) = a so that a n converges towards 1/gi = l/gf'- 
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Table Al. Comparison between the coefficient j n obtained with the Kellogg's method, 
with 0g = r-v/ v(r), and the exact results for the potentials 1(54 (1 -1(56 (1 . 



Potential 


7i 


72 


73 


74 


Exact 


E (£ = 


= 0) 


1.5323 


1.4480 


1.4459 


1.4458 


1.4458 


PE (£ 


= 0) 


0.70463 


0.67718 


0.67669 


0.67668 


0.67668 


SW {£ 


= 0) 


2.4853 


2.4676 


2.4674 


2.4674 


2.4674 


SW {£ 


= 1) 


10.247 


9.8885 


9.8707 


9.8697 


9.8696 


SW {£ 


= 2) 


21.635 


20.317 


20.202 


20.192 


20.191 


SW (£ 


= 3) 


36.630 


33.620 


33.278 


33.227 


33.217 


SW {£ 


= 4) 


55.210 


49.745 


49.000 


48.864 


48.831 


SW (£ 


= 5) 


77.374 


68.669 


67.322 


67.039 


66.954 



Table A2. Comparison between the coefficient /3„ obtained with the Kolomy's 
method, with <fi^ = r^J v (r), and the exact results for the potentials 1)54 (1 - 1)56(1 . 



Potential 


Pi 


ft 


ft 


ft 


Exact 


E (£ = 


= 0) 


1.4674 


1.4465 


1.44582 


1.44580 


1.4458 


PE {£ 


= 0) 


0.68270 


0.67682 


0.67669 


0.67668 


0.67668 


SW {£ 


= 0) 


2.4706 


2.46744 


2.4674 


2.4674 


2.4674 


SW {£ 


= 1) 


10.000 


9.8770 


9.8701 


9.8696 


9.8696 


SW {£ 


= 2) 


20.811 


20.253 


20.198 


20.192 


20.191 


SW {£ 


= 3) 


34.851 


33.439 


33.252 


33.223 


33.217 


SW {£ 


= 4) 


52.1053 


49.372 


48.935 


48.852 


48.831 


SW {£ 


= 5) 


72.5672 


68.022 


67.192 


67.039 


66.954 
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Table A3. Comparison between the upper limit on gc obtain with the variational 
method, the combination of the variational and the Kellogg's method and the exact 
value of <7c . The variational wave function <pe(r) are r exp(— qr) for the E and the PE 
potentials and r p for the SW potential. 



Potential 


<f>£ 




Exact 


E (£ = 


= 0) 


1.44676 


1.44582 


1.44580 


PE (£ 


= 0) 


0.68543 


0.67672 


0.67668 


SW (£ 


= 0) 


2.4747 


2.4674 


2.4674 


SW (£ 


= 1) 


9.9934 


9.8710 


9.8696 


SW (£ 


= 2) 


20.604 


20.201 


20.191 


SW {£ 


= 3) 


34.099 


33.253 


33.217 


SW (£ 


= 4) 


50.357 


48.915 


48.831 


SW {£ 


= 5) 


69.295 


67.117 


66.954 



Table A4. Comparison between the coefficient a n and lu u obtained with the new 
iterative method and the exact results for the potentials (|54|) - (I5GI1 . 



Potential 


a 1 1 






a 4 1 


Exact 


Co> 4 1 


u 3 1 


u 2 1 


-r 1 


E (£ = 


= 0) 


1.3333 


1.4211 


1.4408 


1.4448 


1.4458 


1.4465 


1.4495 


1.4674 


1.6000 


PE {£ 


= 0) 


0.64000 


0.67009 


0.67558 


0.67650 


0.67668 


0.67682 


0.67755 


0.68270 


0.72727 


SW (£ 


= 0) 


2.4000 


2.4590 


2.4664 


2.4673 


2.4674 


2.4674 


2.4677 


2.4706 


2.5000 


SW {£ 


= 1) 


8.5714 


9.4839 


9.7638 


9.8419 


9.8696 


9.8770 


9.9000 


10.000 


10.500 


SW (£ 


= 2) 


15.556 


18.271 


19.452 


19.921 


20.191 


20.252 


20.381 


20.811 


22.500 


SW {£ 


= 3) 


22.909 


28.077 


30.813 


32.149 


33.217 


33.439 


33.801 


34.851 


38.500 


SW (£ 


= 4) 


30.462 


38.506 


43.336 


46.041 


48.831 


49.372 


50.121 


52.105 


58.500 


SW (£ 


= 5) 


38.133 


49.336 


56.671 


61.201 


66.954 


68.022 


69.322 


72.567 


82.500 
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Table A5. Comparison between previously known upper and lower limits J2J, 0. © 
and JHJ on gf* and the exact results for the potentials l|54|l - H56(l . 



Potential 


Eq. © 


Eq. © 


Eq. © 


Eq. © 


Exact 


a 4 1 


u; 4 1 


E (£ = 


= 0) 


1.4383 


1.6755 


1.5442 


1.4467 


1.4458 


1.4448 


1.4465 


PE (£ 


= 0) 


0.67421 


0.76498 


0.86547 


0.67691 


0.67668 


0.67650 


0.67682 


SW (£ 


= 0) 


2.3593 


2.6667 


4.0000 


2.4747 


2.4674 


2.4673 


2.4674 


SW {£ 


= 1) 


9.1220 


11.719 


10.068 


9.9934 


9.8696 


9.8419 


9.8770 


SW {£ 


= 2) 


18.454 


25.413 


20.895 


20.604 


20.191 


19.921 


20.252 


SW {£ 


= 3) 


30.245 


43.570 


35.424 


34.099 


33.217 


32.149 


33.439 


SW (£ 


= 4) 


44.425 


66.089 


53.519 


50.357 


48.831 


46.041 


49.372 


SW (£ 


= 5) 


60.947 


92.909 


75.114 


69.295 


66.954 


61.201 


68.022 



